Energy Spectrum of Spin Fluctuations in Superconducting La2— ajSr^; CuO. 

(0.10 < a; < 0.25) 

Chul-Ho Lee* , Kazuyoshi Yamadat , Yasuo Endoh 
Department of Physics, Tohoku University, Aramaki Aoba, Sendai 980-77, Japan 

Gen Shirane 

Department of Physics, Brookhaven National Laboratory, Upton, NY 11973-5000, USA 
R. J. Birgeneau, M. A. Kastner, M. Greven-'-, Y-J. Kim 
Department of Physics and Center for Materials Science and Enginnering, Massachusetts Institute of 

Technology, Cambridge, MA 02139, USA 

November 12, 1999 



The energy spectrum of incommensurate spin fluctuations in superconducting Sr-doped 
La2Cu04 has been studied by inelastic neutron scattering experiments. An energy gap in the 
spin excitation spectrum is observed in the superconducting state of optimally doped {x=0.15) 
and slightly overdoped (x=0.18) samples. At temperatures well below Tc, the incommensu- 
rate peaks diminish rapidly in intensity with decreasing energy below ~8meV and merge into 
the background below ~3.5meV for x=0.15 and ~4.5meV for x=0.18. For both samples, 
the energy spectrum of the q-integrated spin susceptibility, x"{lo), exhibits an enhancement 
around 7meV, which is caused by a broadening in the momentum width of the incommen- 
surate peaks. The gap-like structure in the energy spectrum and the enhancement of x"i'^) 
survive at Tc. On the other hand, for both underdoped (x=0.10) and heavily overdoped 
(x=0.25) samples there is neither a clear energy gap nor an enhancement of x"{ll)) below Tc. 
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I. Introduction 

The interplay between magnetism and supercon- 
ductivity is one of the most interesting and impor- 
tant issues in the physics of the cupratc high temper- 
ature superconductors [1]. For more than ten years, 
experimental studies on the spin fluctuations in these 
materials have been performed by neutron scatter- 
ing and NMR measurements. NMR first revealed 
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a suppression of the low-energy spin excitations be- 
low what is called the spin gap temperature [2]. In 
the underdoped region, it is thought that above Tc 
a pseudo-gap opens up in the spin fluctuation spec- 
trum. Since the spin-gap state is believed to be re- 
lated to the pairing mechanism, a large number of 
experimental and theoretical studies have focused on 
the origin of the spin gap. 

Unlike NMR, neutron scattering can determine 
both the momentum and the energy dependence 
of the magnetic excitations, providing information 
about the anisotropy of the energy gap of the high- 
Tc cuprates. In fact, neutron scattering has been 
used to observe the spin-gap of YBa2Cu307_y by 
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determining the energy spectrum of the dynamical 
spin susceptibihty around (tTjTt) as a func- 

tion of T [3,4,5]. However, the energy spectrum of 
the spin fluctuations and its doping dependence in 
La2-a;Srj;Cu04 have not yet been determined over a 
wide range of x, most especiahy for x above the op- 
timal doping value. This is in part because it is diffi- 
cult to separate the weak magnetic signals from the 
non-magnetic background originating from phonons 
for example. Further, it is difficult to grow large ho- 
mogeneous single crystals over a wide doping region. 

Several years ago an overview of the low energy 
spectrum of x"i<li'-^) associated with the incommen- 
surate peaks of Lai.85Sro.i5Cu04 was presented [6]. 
Shortly thereafter, a well-defined energy gap of ^ 
3.5meV was observed in the superconducting state 
of the optimally Sr-doped La2Cu04 [7]. The latter 
measurements were limited to low energies ^ 6meV. 
For higher energies, pulsed neutron scattering mea- 
surements were performed on Lai.85Sro.i5Cu04 and 
revealed a broad peak between uj — 40 and 70meV 
[8] in x"(w). The dynamical spin susceptibility was 
integrated over q around the incommensurate peaks. 
This result was recently reconfirmed by more com- 
prehensive experiments and analysis on the same 
crystals [9]. 

However, Hayden et al. [10] reported a somewhat 
sharper peak in the energy spectrum of x"i^) s-t 
around u = 20meV in Lai.ggSro.i4Cu04 by com- 
bining data taken by pulsed neutron scattering with 
that of 3-axis neutron spectroscopy. Furthermore, 
the latter group studied the energy and tempera- 
ture dependence of x"{lli) using the same crystal 
and suggested that they were observing the effects 
of a nearby quantum critical point at T=0 K. They 
reported that the q-width of the spin fluctuations 
monotonically increased with increases in either tem- 
perature or energy [11] although their observed cor- 



relation length at low temperatures was in fact quite 
short in apparent contradiction with the quantum 
critical interpretation. They found no evidence for 
a gap in the spin fluctuation spectrum above Tc, in 
conflict with NMR measurements [12]. 

To reconcile these apparent inconsistencies a more 
systematic neutron scattering study is required. In 
particular, it is important to measure the doping de- 
pendence of the energy spectrum over a wide doping 
region. To accomplish this we have grown large sin- 
gle crystals of superconducting Sr-doped La2Cu04 
over a wide range of Sr concentrations extending 
from the underdoped to the highly overdoped region. 
Although the energy region in the present study is 
still limited to that of conventional low energy 3-axis 
neutron scattering measurements, the doping depen- 
dence provides new insights into the spectrum of the 
spin fluctuations in this system. We describe the ex- 
perimental details in section ||, the results and data 
analysis are presented in section |l|, a discussion is 
given in section ^ and a brief summary is presented 
in section 

II. Experimental Detail 

Single crystals of La2_KSr^Cu04 have been grown 
by the traveling solvent floating zone (TSFZ) 
method using lamp-image furnaces [13,14]. The fur- 
naces have been improved in order to tunc the tem- 
perature gradient near the molten zone. We have 
found that a steep temperature gradient around the 
liquid allows us to grow homogeneous large single 
crystals. As-grown single crystals have been kept 
under oxygen flow at 900°C for 50 hours to remove 
oxygen defects. 

To characterize the Sr content and homogene- 
ity we determined the phase transition tempera- 
ture (Tg) between the high temperature tetrago- 
nal (HTT) and the low temperature orthorhombic 
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(LTO) phases by neutron diffraction. The values 
for Tg of our single crystals [15] agree quite well 
with those previously deduced from powder data 
[16,17,18]. Tg depends primarily on the Sr concen- 
tration and specifically a small amount of oxygen 
non-stoichiometry does not change Tg appreciably 
[19,20]. Therefore, the agreement of the measured 
Tg for our single crystals with the results obtained 
from measurements on powders confirms the Sr con- 
centration of our single crystals. The values of the 
Sr concentration measured by electron probe micro- 
analysis were also consistent with the results from 
the Ts measurements within the instrumental reso- 
lution. 

Superconducting magnetic shielding effects were 
measured by a SQUID magnetometer in a magnetic 
field of 10 Oe. The onset temperature of the su- 
perconducting transition (Tc) in the x = 0.15 single 
crystal, Tc=37.5K, is almost the same as the max- 
imum value found in powder samples [18,21]. For 
X = 0.10,0.18 and 0.25 the values of Tc = 29K, 
36. 5K and 15K, respectively, also agree with those 
of powder samples [15]. Since Tc is very sensitive to 
any deviation from stoichiometry in the oxygen con- 
centration [20,22], this agreement suggests that the 
oxygen content and ordering of our crystals is nearly 
optimal. 

Normal state magnetic susceptibilities were also 
measured using a SQUID magnetometer with a mag- 
netic field of IT. We observed a broad peak at a 
temperature Tmax='^00^ for x=0.18 which depends 
on both the Sr and the oxygen concentration [22]. 
Tmox of the single crystal is in good agreement with 
that of a powder sample with the same Sr concen- 
tration [22, 23, 24]. Lattice parameters of pulverized 
single crystals measured by x-ray diffraction at room 
temperature are also consistent with previous data 
[15]. 



Inelastic neutron scattering measurements 
were performed with the triple-axis spectrometer, 
TOPAN at JRR-3M in the Tokai EstabUshment of 
JAERI for crystals with x=0.15, 0.18, and 0.25 and 
with H7 at the HFBR at Brookhaven National Lab- 
oratory for crystals with x = 0.10 and 0.15. The 
incident (final) neutron energy was typically fixed 
at Ei (E/) = 14.75meV for TOPAN and Ej (Ef) = 
14.7meV for H7. A typical sequence of the horizon- 
tal collimators was 40'-100'-S-60'-B for TOPAN and 
40'-40'-S-80'-80' for H7 where S denotes the sam- 
ple position. Pyrolytic graphite crystals were used 
both as monochromator and analyzer. A pyrolytic 
graphite filter was used to reduce the intensity of 
higher order neutrons. Moreover, in TOPAN, a sap- 
phire crystal was inserted to reduce significantly the 
flux of high energy neutrons. The single crystals 
were mounted in an Al container filled with He gas 
as a heat exchanger. For cooling, a closed cycle He 
refrigerator was used. 

III. Results 

The four panels in Fig. 1 show representative q- 
spectra of the magnetic fluctuations measured at and 
below Tc. Constant-energy scans were performed 
through the two peaks at (tt, 7r(l ± S)) as illustrated 
in the inset of panel (c). For x = 0.18, due to some 
non-magnetic contamination in the scan, we present 
only one of the peaks, at (tt, 7r(l — 6)). For T ~ 
Tc well-defined peaks are observed at incommensu- 
rate positions for all samples. Below Tc, the peaks 
at 3 meV and 3.5 meV vanish into the background 
for x = 0.15 and 0.18 respectively. As is shown by 
the energy spectra discussed below, this dramatic 
change occurs due to the opening of an energy gap 
in the spin excitation spectrum at the same tem- 
perature as that for the onset of superconductivity. 
In contrast, for x = 0.10 and 0.25 substantial peak 
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FIG. 1. Constant energy-scan across {iv,n) as indicated 
in the inset of panel (c) for samples with (a) x=0.10 (b) 
x=0.15 (c) x=0.18 and (d) x=0.25 below and at T^. All 
solid lines and dashed lines for x = 0.10 and 0.25 are the 
results of least square fits to the data points using two 
Lorentzian functions convoluted with the instrumental 
resolution. 



intensities at low energies are observed even at tem- 
peratures well below Tc 

We present the results of the same q-scan as that 
shown in Fig. 1 at several different energies for 
x=0.10 and 0.15 in Figs. 2 and 3 respectively. Be- 
low Tc, the peak intensities for x=0.10 are compa- 
rable at all energies between 2 and 10 meV. On the 
other hand, for x=0.15 the intensities dramatically 
decrease with decreasing energy and no peak remains 
outside of the errors at 3 meV. Although our collabo- 
ration has reported the observation of an energy gap 
for x=0.15 previously [7] the change in peak inten- 
sities below ~7 meV appears to be more dramatic 
in these experiments. A similar energy dependence 
of the incommensurate peak intensity is found for 
x=0.18 as shown in Fig. 4(c). 

The q-spectrum of the incommensurate peaks was 
fitted using the following dynamical structure factor 
S(q,ti;) convoluted with the instrumental resolution 
function: 

/ 



S{q,uj) = 



1 — exp{- 



kbTI 



with 



5=1.4 



\q - qs\^ + nl 



(1) 



(2) 



where qg,K^,kB and are the peak positions of 
the magnetic peaks around (7r,7r), the q-width at a 
given w, the Boltzman constant and an overall scale 
factor, respectively, assumed to be isotropic in 

the Cu02 plane. 

The imaginary part of the q-integrated dynami- 
cal susceptibility averaged over the Brillouin zone is 
defined as: 



X"{cj)/u;^ / x"{q,u;)dq/{u;- / dq) 



(3) 

Note that in this case, the q-intcgration is carried 
out for the four incommensurate peaks around (tt, tt) 
over the Brillouin zone. 
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FIG. 2. Q-spectra at different energies obtained by the 
same scan as in Fig. 1 for x=0.10 below Tc 
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FIG. 3. Q-spectra at different energies obtained by the 
same scan as in Fig. 1 for x=0.15 below Tc. 
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FIG. 4. Energy dependence of the intensity at the in- 
commensurate peak position below Tc for a;=0. 10(a), 
0.15(b), 0.18(c) and 0.25(d). Panels (e) - (h) depict the 
energy dependence of x"{^)/^ calculated by using the 
fitted results of the q-spectra. 

In Fig. 4, we show the energy-dependence of 
the peak intensities at the incommensurate positions 
for X = 0.10,0.15,0.18 and 0.25 together with the 
resolution-corrected which is obtained by 

fitting the q-spectrum and integrating over q for the 
four incommensurate peaks around (tt, tt). The abso- 
lute value of x"(ijj)/ti> for each sample was calculated 
by using phonon intensities for the same sample as 
described in the Appendix. We present x"(a))/a; in- 
stead of x"(cli) to reveal any gap structure in the 
spectrum because the latter quantity may go to zero 
as the energy approaches zero even if there is no 
gap. Fig. 4 again exhibits two types of energy depen- 
dencies for the incommensurate peaks for the four 
samples. There exists a well-defined energy gap for 
X ^ 0.15 and 0.18, accompanied by an enhancement 
of x"{lo) /uj around 6 meV to 7 meV while no gap-like 
structure is seen for x — 0.10 and 0.25. 

The energy dependence of the resolution-corrected 
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FIG. 5. Energy dependence of the q-width of an incom- 
mensurate peak below Tc with the instrumental resolu- 
tion deconvolved. The q-width in the low energy region 
for x=0.15 and 0.18 cannot be obtained due to the open- 
ing of the energy gap below Tc. 

q-width of the incommensurate peak below Tc is 
shown in Fig. 5. We only show the data above 5 
and 6 meV for x = 0.15 and 0.18, respectively be- 
cause of the energy gap. We note that the q-width 
has a rather different lu dependence depending on 
whether or not there is a gap. For the gapped sam- 
ples, above the gap energy the q-width is small and 
it increases with decreasing energy as the gap energy 
is approached. Furthermore, as will be shown later 
in Fig. 8, for T the q-width appears to exhibit a 
peak for excitation energies near the gap energy. On 
the other hand, for the gapless samples the q-width 
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FIG. 6. Energy dependence of x"{^) for x=0. 15(a) and 
0.18(b) around Tc. The solid lines are the results of 
a least squares fit using the gap-function, Eq. (4), con- 
voluted with the instrumental resolution function. The 
broken line in the upper panel is a smoothed interpolta- 
tion of the data for x — 0.14: reported in ref. [10]. 

is constant with energy at least up to 12 meV. 

Next, we show x"{^) for temperatures T^Tc for 
La2_xSrxCu04 crystals with x=0.15 and 0.18 in Fig. 
6. Although the peak in the intensity appears at 
slightly lower energies compared to the data be- 
low Tc, a gap- like structure and an enhancement in 
x"(w) still remains. In the upper panel of Fig. 6, we 
show using a dashed line the data obtained by Hay- 
den et al. for the same system with x=0.14 [10]. No 
peak in x" {^) is observed in their data. As shown in 
Fig. 7, for the gapless sample of La2-xSrxCu04 with 
x=0.10, both the energy and temperature variation 
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FIG. 7. Energy spectra of x" for x=0.10 below and 
around Tc. 

of is weak, although we see a slight suppres- 

sion of x"{^) around 4 to 6 meV at T=9K. 

In Fig. 8 we show the energy variation of the 
q-width of the incommensurate peak for T^^T^ for 
x=0.15 and in the inset for x=0.18. In the figure, we 
have plotted results from several different crystals. 
We observe a peak in the q-width at around the same 
energy as that at which the enhancement of x"('-^) 
is observed. In contrast, the energy variation of the 
q-width for the sample with x=0.14 of Aeppli et al. 
[11] shows a linear energy dependence. It should be 
noted that the peak width for energies of order 8 
to 10 meV is substantially narrower below Tc than 
that at Tc for the gapped samples, particularly for 
x=0.18. 

IV. Discussion 

The present study has revealed two types of en- 
ergy spectra for the spin fluctuations in Sr-doped su- 
perconducting La2Cu04. For x — 0.15 and 0.18, a 
well-defined energy gap appears below Tc, whereas 
for the underdoped sample, x — 0.10, and highly 
overdoped sample, x=0.25, no gap is observed even 
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FIG. 8. Energy dependence of the q- width of an in- 
commensurate peak around Tc deconvoluted with the 
instrumental resolution for x=0.15 and 0.18(inset). The 
broken line is drawn to represent the data for x=0.14 
reported in ref. [11]. 

though they exhibit bulk superconductivity with Tc 
=29K for a; = 0.10 and T^^ISK for x=0.25. 

In order to discuss the resuhs shown in Figure 
4 quantitatively, we introduce a phenomenological 
dynamical spin susceptibility of the form: 



B- 
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7^ + uj^ 



Re 



Lo — iV 



(4) 



This corresponds to Lorentzian spin fluctuations 
in the normal state multiplied by a gap-function in 
the superconducting state. The gap function in- 
cludes the gap-energy A together with a broadening 
F; 7 is the inverse of lifetime of the spin fluctua- 
tions. The parameters B, A and F are assumed to 
be constant in the q-region of the incommensurate 
peaks. It should be noted that in Eq. (4), 
is a simple Lorentzian in energy for A = 0. On the 
other hand, in the superconducting state, Eq. (4) 
describes the gap-structure with a non-zero value of 
A. As required, the above dynamical spin suscep- 
tibility satisfies detailed balance for the dynamical 
structure factor S(q, lo). 



The observed energy spectra were fitted using Eq. 
(4) convoluted with the instrumental energy resolu- 
tion. The value of A so-obtained is 6 ± 0.3 meV 
for both x=0.15 and 0.18 while F = 0.6 ± 0.3 meV 
for x=0.15 and = ± 0.2 mcV for x=0.18. Quite 
recently. Lake et al. reported a similar value for 
the energy-gap for a sample of La2-xSrxCuOi with 
x=0.16 [25]. We note that A is significantly larger 
than the value previously quoted by our collabora- 
tion which was around 3.5 meV for x=0.15 [7]. How- 
ever, this latter value for the energy-gap was simply 
the energy below which the magnetic intensity dis- 
appeared entirely below Tc. This energy is naturally 
less than our fitted value simply due to the non-zero 
broadening of the gap-structure for x=0.15 together 
with the effects of instrumental resolution. We note 
that for samples with x— 0.14 and 0.16, Lake et al. 
[25] find F = 1.2 meV and 0.1 meV respectively al- 
beit using a somewhat different functional form than 
Eq. (4). Thus, F decreases smoothly with increasing 
X reaching within the errors for x=:0.18 while the 
gap energy A appears to be independent of x for x 
between 0.14 and 0.18. 

In Ref. 6 the gap-energy, A™'^ of the spin fluc- 
tuations at the incommensurate peak position was 
converted into the value at (7r,0) by assuming a d{x^- 
y2)-type q-dependence. The previous value A™'^ 
'--^3.5 meV implied a full gap of ^10 meV which is 
consistent with the full superconducting gap of ^^8 
meV inferred from photoemission spectroscopy [26] . 
However, the new value of A™^ determined here im- 
plies a full gap of ^^18 meV which corresponds to 
A = b.QkTc which seems too high compared with 
the BCS weak coupling value of 1.77 kTc. This fact, 
therefore, suggests that A™^ may not be simply re- 
lated to the superconducting gap. It is possible that 
A™^ instead relates to the so-called pseudo-gap or 
the larger energy gaps reported by photoemission 
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and tunneling spectroscopies, 

The gapped samples exhibit an anomalous broad- 
ening in the q-width of the incommensurate peaks 
for energies around the gap value of AJ^'^. Qualita- 
tively, we can interpret the peak-broadening as aris- 
ing from the energy independent dispersion surface 
of the quasiparticle excitations around the energy- 
gap. In other words, if the magnetic scattering oc- 
curs through interband excitations of quasiparticles, 
the momentum change can be widely distributed for 
excitation energies close to the gap. As shown in Fig. 
8, a broadening in q of the peak width for energies 
of order A is observed at Tc which again indicates 
the robustness of the gap structure of x"(w) for the 
gapped samples. 

We discuss next the possible reasons for the ab- 
sence of an energy-gap in the gapless samples with 
x=0.10 and x=0.25. At present, it is widely specu- 
lated without any direct experimental evidence that 
the disorder introduced by dopant substitution at 
the La-sites degrades the lifetime of the quasiparti- 
cles and smears out the gap structure. However, this 
is contradicted by the fact that no gap is observed 
in either stage-6 or stage-4 La2Cu04+y in spite of 
the fact that structural disorder is minimal in these 
systems [27,28]. Further, neutron scattering experi- 
ments on the Yl-2-3 system reveal the robustness 
of the magnetic energy-gap against Zn-impurities 
[29,30]. In the Zn-doped Yl-2-3 system, although 
low energy spin fluctuations appear at low tempera- 
tures, the gap-like or pseudo-gap structure appears 
at higher temperatures. 

For the x=0.10 sample, the low energy spin exci- 
tations are simply explained as originating from the 
spin wave excitations concomitant with the elastic 
magnetic order which occurs in the superconduct- 
ing state. Such order is observed for all x between 
0.02 and 0.135 in La2Cu04. However, a different ex- 



planation is required for the highly overdoped gap- 
less sample, x=0.25, because in that case there is 
no evidence for magnetic order coexisting with the 
superconductivity. Here, we simply point out that 
there is a dramatic degradation of the spacial coher- 
ence of the spin correlations in the overdoped region 
particularly near the upper critical doping value for 
the superconductivity as pointed out by Yamada et 
al. [15]. In this doping region, the spatial coherence 
length for the low energy spin fluctuations becomes 
comparable to or shorter than the superconducting 
coherence length. Therefore, the energy gap in this 
region is probably vitiated by the degradation of the 
lifetime of quasiparticles due to the short ranged spin 
fluctuations. 

V. Summary 

We have studied the energy dependence of the spin 
fluctuations for superconducting La2-xSrxCu04 
(a;=0.10, 0.15, 0.18, 0.25) using neutron inelastic 
scattering. Below Tc a well-deflned energy gap in 
the incommensurate spin fluctuations is observed in 
the superconducting state of the optimally doped 
(a; = 0.15) and slightly over doped {x = 0.18) sam- 
ples. An enhancement in x"{^) caused by the peak- 
broadening is also observed at ~6 meV which re- 
mains even at Tc. For the underdoped sample with 
X = 0.10 and the highly overdoped sample with 
X = 0.25, no clear gap is ovserved even though these 
samples show bulk superconductivity. 
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Appendix 

We have converted the observed magnetic inten- 
sity into the absolute value of x"(a') by comparing 
with the intensity of phonon scattering. The phonon 
intensity is given by 



(5) 



where a is the detector efficiency and I„ is the 
energy-integrated phonon intensity measured by 
constant-Q scans. The energy integrated phonon 
scattering cross section is given by: 



\dn) ^ M 



e < n+ 1 > 
(6) 



where C=2.09 barns, F is the nuclear structure fac- 
tor, M is the total atomic mass in a unit cell, (3 is 
the angle between Q and the phonon polarization 
direction, N is the number of unit cells, e~^^ is the 
Debye- Waller factor and < n+ 1 > is the population 
factor. We have measured the integrated inensity 
of the acoustic phonons by constant-Q scans at (2, 
-0.13, 0) for X = 0.10 and (2, -0.17, 0) for x = 0.15, 
0.18, and 0.25. The magnetic scattering intensity is 
measured by constant-w scans. To scale the mag- 
netic scattering intensity with the phonon intensity 
measured in a different scan mode, is divided by 
the sound velocity. The absolute value of x"{q,uj) 
is estimated using the following magenetic cross sec- 
tion per formula unit for a Heisenberg system with- 
out Neel ordering: 



\dflduj J 



„2 



J2^l-kl)S'^{Q,u;) (7) 



5«(g,w) = -^i^ <n+l> Jmx"(Q,u;) 



(8) 



Here Tq = 0.291 barns, f{Q) is magnetic form factor 
and k = Q/\Q\. 
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